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ABSTRACT 


A theoretical  study  of  ice  crystal/ shock  layer  interaction  has  been 
carried  out  in  order  to  assess  the  relative  kinetic  energy  of  an  ice  crystal 
upon  impact  with  a hypersonic  vehicle.  A dynamic  fracture  criterion  has 
been  developed  which  predicts  the  size  of  the  ice  crystal  fragments  formed 
at  the  vehicle  bow  shock.  An  analysis  of  the  ice  crystal  breakup  predicts 
that  the  outer  fragments  continue  to  shatter  and  prevent  the  shock  layer 
gas  from  penetrating  the  fragment  cloud.  The  fragment  cloud  is  then 
treated  as  a deformable  fluid  drop,  and  the  reduced  impact  velocity  due 
to  lateral  spreading  is  determined. 
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I.  INTRODUCTION 


Asa  supersonic  or  hypersonic  vehicle  passes  through  a real  atmos- 
phere. it  may  encounter  relatively  large  clouds  of  dust,  rain  or  ice.  A 
great  deal  of  attention  has  previously  been  given  to  the  problems  of  dust^ 
and  rain^"*1^,  but  not  until  recently  has  significant  effort  been  devoted  to 
the  interaction  of  ice  crystals  with  the  shock  layer.  has  studied  the 

melting  and  vaporization  of  an  ire  crystal  in  a hypersonic  shock  layer  and 
demonstrated  a significant  reduction  in  particle  size  and  impact  velocity 
for  ice  crystals  initially  of  the  order  of  10M  diameter.  Larger  crystals 

were  ralatively  unaffected  by  melting  and  vaporization.  In  the  same  report, 
(12) 

Tea  re  demonstrated  that  ice  crystals  may  shatter  at  the  vehicle  bow- 
shock.  This  introduces  the  possibility  that  the  fractured  ice  crystal  or 
"fragment  cloud"  may  deform  as  it  passes  through  the  vehicle  shock  layer. 
This  deformation  will  be  lateral  to  the  relative  gas  velocity  and  will  tend  to 
increase  the  drag  area  of  the  fragment  cloud.  Hence,  the  velocity  of  the 
fragment  cloud  will  decrease  more  rapidly  than  that  of  a rigid  particle,  and 
the  velocity  with  which  the  fragment  cloud  impacts  the  vehicle  may  be  signi- 
ficantly lower  than  that  with  which  a "rigid"  ice  crystal  would  impact  the 
vehicle. 

To  study  this  phenomenon  we  must  first  develop  a criterion  for  the 
dynamic  fracture  of  the  ice  crystal  at  the  vehicle  shock.  This  is  carried  out  in 
section  II.  The  problem  of  ice  crystal  breakup  is  addressed  in  section  HI, 
and  in  section  IV  we  examine  the  dynamics  of  the  fragment  cloud  and  demon- 
strate the  reduced  impact  velocity  due  to  the  lateral  deformation  of  the  cloud. 
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The  dynamic  shattering  criterion  used  in  section  II  is  subject  to  the  experi- 
mental verification.  However,  the  dynamics  of  the  fragment  cloud  are  insen- 
sitive to  the  actual  shattering  process  and  the  results  of  section  IV  are  only 
weakly  dependent  on  the  dynamic  shattering  criterion. 


II.  INITIAL  FRACTURE  OF  THE  ICE  CRYSTAL 


A)  Conditions  for  Fracture 

A schematic  of  an  i : crystal  entering  a hypersonic  shock  layer  is 
given  in  Fig.  1.  The  ice  crystal  enters  the  vehicle  shock  layer  with 
velocity  U , finds  itself  moving  with  supersonic  velocity  relative  to  the 

49 

shock  layer  gas  and  thereby  generates  its  own  bow  shock  or  crystal  shock. 
Only  the  normal  component  of  the  shock  Mach  number  is  important  for  the 
pressure  loading  on  the  crystal.  The  pressure  loading  as  a function  of  normal 
Mach  number  can  be  calculated  with  the  aid  of  Fig.  2.  The  subscript  x 
represents  conditions  upstream  of  the  vehicle  shock  and  subscript  y denotes 
conditions  in  the  vehicle  shock  layer.  The  Mach  number  of  the  ice  crystal 
relative  to  the  shock  layer  gas  is 


which  may  be  determined  from  the  normal  shock  tables.  The  pressure  Pg 

is  the  total  pressure  behind  a shock  of  Mach  number  Mj  moving  into  a gas 

of  pressure  Py.  Thus,  P8  may  also  be  obtained  from  the  normal  shock 

13 

relations.  The  static  yield  criterion  indicates  that  fracture  will  occur  in  the 
vehicle  shock  layer  if  the  radius  of  Mohr's  circle,  (Ps-Py)/2,  is  greater  than 
Ss,  the  yield  stress  of  ice  in  shear. 


If  the  speed  of  sound  in  the  ice,  a0,  is  greater  than  Vx,  the  ice  crystal  will 
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ICE  CRYSTAL 
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VEHICLE  SHOCK 


RYSTAL  SHOCK 


VEHICLE 


% 1 Ice  Crystal-Shock  Wave  interaction 


shatter  upstream  of  the  vehicle  shock  if 

P -P  > 2 S . 

S X — 8 

This  fracture  criterion  is  illustrated  in  Fig,  3 and  demonstrates  that  most 
hypersonic  vehicles  will  shatter  ice  crystals  through  the  bow  ,hock  wave  - 
ice  crystal  interaction  process  outlined  above.  However,  this  conclusion  is 
based  on  the  application  of  a static  yield  criterion.  We  must  proceed  to 
develop  a model  for  the  dynamic  fracture  of  a brittle  material  under  impulsive 
loading. 

B)  Dynamic  Fracture  Criterion 

Ice  crystals  passing  the  bow  shock  wave  of  a hypersonic  vehicle  are 
subject  to  very  short  duration  impulsive  loads,  and  a consistent  dynamic  yield 
criterion  must  be  established  in  order  to  analyze  the  shattering  problem. 

First,  we  define  three  different  times  relevant  to  dynamic  loads.  Tp)  is  the 
duration  of  the  applied  load,  Ty  is  the  transient  time  for  application  of  the 
load  (Ty  <.  Tq)  and  Tj^  is  the  response  time  or  the  time  required  for  a 
specimen  of  length  L to  know  a load  has  been  applied.  Clearly,  = L/a.^, 
where  aQ  is  the  speed  of  sound  in  the  specimen. 

By  definition,  static  loading  occurs  when  Ty  » and  the  static  yield 
stress  (Oc)  is  the  appropriate  measure  of  the  maximum  load  that  the  specimen 
can  tolerate  without  plastic  deformation  or  failure.  Clearly,  when  Ty  « t^, 
a dynamic  yield  stress  (o  j.)  is  a necessary  measure  of  failure.  When 
= O { i _)  we  assert  ’■hat  o ^ = O (a  ). 

1 K D 8 


-6- 


UPSTREAM  PRESSURE  (ATM) 


fig.  3 Conditions  for  Fracture 


One  example  of  a situation  where  is  a monochromatic 

stress  wave.  Let  Xx  be  the  longitudinal  stress  (tension  or  compression) 

in  a stress  wave 


X = A.  sin  (nka  t)  cos  (TTkx). 

X x ® 

Therefore,  failure  occurs 
when  - 0 (og).  An  alternate  way  of  looking  at  the  stress  wave  is  in 


For  an  element  of  length  1/k,  tr  = ^ 


terms  of  the  particle  velocities,  U . The  momentum  equation 


P 


X 

x 


yields 


a Av 

U = — ~ — cos  (ftka  t)  sin  (nkx) 
x Jt-  o 


where  E is  young's  modulus  relating  stress  (a)  to  strain  (e) 


(o  = 


Ee  and 


*0=  rfETp)-  If  w«  interpret  Ux  as  a series  of  elements  of  length  1/k 

(14) 

colliding  with  velocity  aQA^/E,  then  impact  theory  indicates  that  failure 
occurs  when  the  impact  velocity  exceeds  C (e  a ).  Since  o = Ee  , this 

' SO  8 S 


criteria  reduces  to  A,  = 0 (o  ) 

k s 

Having  developed  a rational  argument  for  applying  a static  yield 
criterion  to  the  monochromatic  stress  wave,  we  now  make  the  following 
bold  hypothesis  which,  while  plausible,  c*n  be  justified  only  if  the  subsequent 
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conclusions  agree  with  experiment. 


HYPOTHESIS:  Dynamic  Yield  St  re  a * may  be  obtained  by  Fourier 

expending  the  applied  load  end  applying  the  monochromatic 

yield  criterion  to  each  wave.  Failure  occur*  when  A,  = Q 0 

k • 

where  Q is  O (1)  and  ie  to  be  determined  by  applying  this 
technique  to  an  impulsively  loaded  specimen. 

The  primary  argument  agtinat  this  hypothesis  is  that  the  single  monochro- 
matic wave  may  in  no  way  represent  the  total  stress  which  is  the  sum  of  all 
the  waves.  While  this  is  true,  one  could  assert  that  a specimen  of  length  L 
would  be  fractured  by  its  primary  mode  and  that  a small  segment  (length  4) 
of  the  specimen  would  be  shattered  by  the  primary  mode  associated  with  the 
length  4,  even  though  the  primary  mode  of  the  small  segment  would  be  a 
higher  order  mode  of  the  larger  specimen.  Hence,  it  is  within  this  context 
that  we  apply  the  monochromatic  yield  criterion  to  each  Fourier  component 
of  a dynamic  load.  While  we  cannot  theoretically  confirm  the  basic  hypothesis, 
we  proceed  with  the  intention  of  verifying  it  through  empirical  confirmation 
of  the  conclusions  which  are  subsequent  to  the  basic  hypothesis. 

C)  One -Dimensional  Fracture 

Let  us  first  examine  the  present  theory  by  applying  it  to  the  problem 

of  a one -dimensional  slab  of  thickness  b that  is  impulsively  loaded  on  one 

side  with  pressure  Pg  for  time  The  load  is  impulsively  removed  st  the 

end  of  that  time  interval  (t_  « b/a  ).  A schematic  of  the  loading  and  un- 

D o 

loading  waves  in  (x,  t)  coordinates  is  given  in  Fig.  4.  If  we  stand  at  some 
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3b /a 


Fig.  4 Wave  Diagram  of  I mpulsive 
Short  Term  Load 


position  x between  *oTp/2  end  b-a^p/2,  the  time  dependent  force  on  en 
element  is  shown  in  Fig.  5.  The  Fourier  expansion  of  is  given  by 


The  dynamic  fracture  criterion  indicates  that  the  primary  mode  {n  = 1)  will 
fracture  the  specimen  if 


> Q O 

8 


or 


Hence, the  dynamic  yield  stress  is  given  by 


and  illustrates  that  Op  is  inversely  proportional  to  and  indeed  is  in  excess 
of  o#  provided  that  Q is  of  order  unity. 

To  determine  the  value  of  Q,  let  us  again  take  a one -dimensional  slab 
and  impulsively  load  it  in  tension  (T)  for  all  time.  The  wave  diagram  is 
illustrated  in  Fig.  6 and  the  pres  sure -time  history  is  shown  in  Fig.  7.  The 
corresponding  Fourier  series  representation  of  the  dynamic  load  is  given 

i 

I 

by 


-11- 


Fig.  5 Pressure  "Trace”  For  Impulsive 
Short  Term  Load 


Fig.  6 Wave  Diagram  For  Impulsive  Load 
in  Tension 
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X - T +'N"'1  — (co*  ntT  - 

x 2L#  nTT  \ 

n 

An  impulsively  applied  load  creates  twice  the  strain  as  a slowly  applied 
load.  Thus,  the  impulsive  load  will  fracture  the  specimen  when  T is 

only  0(/2.  Equating  the  amplitude  of  the  n = 1 mode  of  the  above  expansion 
to  QOg*  we  determine 


Note  that  we  have  obtained  Q by  comparing  the  amplitude  of  standing  waves. 
Any  subsequent  use  of  the  value  of  Q must  be  made  with  standing  w&ves.  A 
completely  consistent  approach  could  be  developed  using  traveling  waves  but 
the  values  of  Q would  differ  by  a factor  of  2. 


Having  developed  a criterion  for  fracture  of  a dynamically  loaded 
specimen,  let  us  now  examine  shattering.  Consider  a one -dimensional  slab 
impulsively  loaded  on  one  side  and  held  constant  for  all  time.  This  is  the 
one -dimensional  analogue  of  the  ice  crystal  loading  at  the  bow  shock  of  a 
hypersonic  vehicle.  The  Fourier  expension  of  the  dynamic  load  is  expressed 
as 


Suppose  P »o  . Then,  not  only  can  the  n = 1 mode  fracture  the  specimen, 
s * 

but  \v.he  n = 2 mode  can  break  those  fragments  and  so  on  up  to  the  critical 
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n 


▼aloe  of  n (bq)  where  the  nth  mode  is  just  strong  enough  to  fracture  the 
specimen.  This  value  of  n is  determined  fay 


= Qo 

n tr  s 

c 


n ~ — 

c <j 


and  the  corresponding  sire  of  the  fragments  (1)  of  the  n mode  is  give^  by 

c 


■ft)-  - 


These  are  the  largest  fragments  that  can  survive  the  applied  pressure  P . 
Smaller  fragments  can  exist  due  to  the  shattering  by  other  modes.  That  is. 
the  n^th  mode  fractures  the  specimen  ar  values  of  x/b  equal  to 


x JL  _3_  _5_  _7_ 

b = 2n  ’ 2 n ’ 2n  ’ 2 n 
c c c s 


2 n - 1 
c 

2 n • 
c 


The  n^-i  mode  fractures  th  t specimen  at 


b ' 2 <n  -1)’  2 (n  * 1 ) ’ 2 (n  -1)’  2 (n  -1) 

c c c c 


2 (n  - 1)  - 1 
c 

2 (n  - 1) 
c 


The  largest  fragment  is  then 


2 n 2 n 
c c 


wh«r«u  the  imalleit  is 


&x  _ 1 1 

b 2 (n  .1)  ' 2n 


2 n 


Therefore,  for  one -dimensional  shattering,  we  have  determined  that  the 
fragment  sixes  range  from  the  smallest 


mm 


i 

rip) 


to  the  largest,  i 


max 


max 


(2) 


Since  the  ice  crystal  erosion  problem  is  dominated  by  the  large  crystals, 
we  shall  study  the  two-dimensional  and  three-dimensional  shattering 
problems  with  the  understanding  that  the  characteristic  fragment  size  refers 
to  the  largest  fragments,  and  that  smaller  fragments  can  exist, 

D.  High  F rcqucncy  Limitations 

The  results  of  section  C)  places  no  upper  limit  on  the  frequency  of  the 
fracturing  wave  or,  equivalently,  place  no  lower  limit  on  the  size  of  the 
shattered  fragments.  It  is  well  known^*^  that  the  minimum  ™ave  length 
that  can  exist  in  a solid  is  of  the  order  of  the  lattice  spacing  (s)  and  the 
maximum  frequency  is  the  photon  frequency  u> 


il 


uj  ■■  a ft 

P 0 
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\ 


mm.  •!•**■  *: 


where  »o  is  the  speed  of  sound  in  the  specimen*  These  concepts  may  be  in- 
eluded  is  the  shattering  analysis  by  adding  some  "inertia"  to  the  stress  tensor. 


Assuming 


o = E(e+M) 


where  e is  ~~  and  both  £ and  0 are  constants*  the  momentum  equation 
Ot 


becomes 


2 2 3 

„ & u ao  _ b u . __  b u 

bt2  **  J*2  bthx2 


Solving  the  initial  value  problem. 


U (x, t)  = sin 


for  X real  and  m complex  we  see  that  propagation  cannot  exist  above  the  fre- 
quency 2/0  or  below  the  wavelength  Tta^B.  Hence,  it  follows  that 


Solving  the  boundary  value  problem. 


-x/6  i (uu  t - 2 ttx />.  ) 


U(x,t)  = e 


for  uu,  A.  s-nd  the  decay  depth  6^  all  real,  we  obtain 


2fi  2 

tt)  P 8 JU 


where  b ^ approaches  s as  tu  approaches  x>  . 


For  an  applied  load  of  arbitrary  frequency  uu,  caution  must  be  takr  n 


to  insure  that  > L before  using  the  proposed  dynamic  yield  criterion  since 
it  omit'  skin  depth  effects.  For  an  impulsively  loaded  specimen  of  length  L, 
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tht  frequency  corresponding  to  the  n'th  mode  is 


n a 

o 


Validity  of  the  shattering  analysis  requires 

6 > L 

e 

or  L > n2s 


or  l > ns 

max 

and  L . > s 

nmn 

Thus,  the  ice  crystal  shattering  analysis  is  valid  down  to  fragment  sizes  of 
the  order  of  s.  High  frequency  effects  do  not  influence  the  shattering  by  an 
impulsively  applied  load  but  may  become  important  for  high  frequency  loads 
and  large  test  specimens. 

E.  Two-Dimensional  Fracture 

The  one -dimensional  shattering  analysis  is  readily  extended  to  two 
dimensions  by  expanding  the  applied  load  into  a double  Fourier  series. 
Suppose  we  have  a two-dimensional  ice  crystal  with  pressure  Pg  impulsively 
applied  at  one  end.  The  applied  pressure  may  be  expanded  as  a square  wave 
over  the  lateral  dimension  (a)  of  the  ice  crystal 


The  one -dimensional  plane  wave  solution  was  given  by 


whe  re  a wag 
o 


^E7p  and  P#  (y)  wa»  constant.  Thus,  the  plane  wave  represen. 
tation  of  an  end  loaded  two-dimensional  ice  crystal  is  given  by 


m n 


-4  Pj  (1  -cos  m tt) 


am 


m 


mntT 


\ ( \ /nrra  t^ 

a J9in[  b jCO#y  b y 


or,  in  standing  wave  form 


• ZZ 

m n 


-2  P (l-cosmrr)  /nira  t 


•cos 


mrr 


cos 


lE 


ttx  mrry  ( n 

- — - u > - cos  < — 

b . ( \ 


rrx  m 
b a 


H] 


(3) 


Similarly,  we  take  a two-dimensional  specimen  ?nd  impulsively  load  it  in 
tension  T for  all  time. 


- \ ' \ ' 2 T (1  - cos  mTT)  (cos  nn-1) 


x L-t 

m n 


cos 


nrra  t 
o 


mn  tt 


1 

l 

- 

• 

cos  | 

1 nrrx 

m TTy  ) 

[ c0fifnTTX  + mTTy  1 

( b 

a 

i » 

_ 

/ 

\ ' 

Just  as  in  the  one -dimensional  case,  failure  occurs  when  the  tension  T 
exceeds  one  half  of  its  static  limit.  Therefore,  the  coefficient  relating  the 
wave  amplitude  to  the  yield  stress  is 


D - W * 


14) 


The  value  of  *8  no*  rest:r'cte^  to  * particular  failure  mode.  The 
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maxim  am  tension  T may  be  limited  by  either  shear  or  tension.  In  obtaining 

Eq.  (4),  it  w*s  stated  only  that  the  impulsive  load  is  twice  as  strong  as  the 

2 

static  load.  The  primary  mode,  whose  amplitude  is  8 T/n  , will  break  when 
T is  "1/2".  Hence,  Q2D  is  4/tT2. 

For  an  axial  load  such  as  that  applied  to  die  ice  crystal,  Mohr's  circle 
predicts  that  the  maximum  shear  (S)  occurs  at 


and  is  of  strength 


y = + x v constant 


Examining  Eq. 
along 


(3),  we  see  that  the  strength  of  the  (m,  n)  mode  is  a maximum 


y 


+ constant 


Therefore,  the  natural  mode  selection  that  would  maximize  the  shear  stress 
on  the  wave  is 

a _ m 
b ~ "n 

Under  the  constraint  of  b a a,  hence  n a m,  the  shear  stress  on  the  two- 
dimensional  ice  crystal  may  be  expressed  as 
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The  dynamic  fracture  criterion  predicts  failure  for  the  critical  value  of  m. 


denoted  by  m . Equating  the  wave  amplitude  of  S to  S , 

2 P a 2D  * 

TT  = s.Q2Dtor(b^*>- 
m It  b 
c 

the  value  of  rr^  becomes 


m 


.=V(5T5)’  ft-*  0- 


(5) 


Note:  for  m < 1,  we  must  use  the  one -dimensional  res  ults  of  Section  C, 

otherwise  the  lateral  dimension  of  the  fragment  is  greater  than  a.  The 

results  are  continuous  if  it  is  assumed  that  the  maximum  axial  stress  is 

limited  by  the  shear  mode  (Og  = 2 S^).  The  characteristic  fragment  size 

is  a/m  or 
c 


Hence,  the  two-dimensional  shattering  analysis  is  appropriate  only  if  the 
fragment  size  is  lest  than  the  small  dimension  of  the  initial  crystal.  Other- 
wise, shattering  is  one -dimensional  along  the  length  of  the  specimen. 

F . Three  -Dimensional  Shattering 

Suppose  that  we  now  have  a three-dimensional  crystal  of  length  b 
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*nd  base  area  a (b  a a).  Expanding  the  applied  force  in  a triple  Fourier 
aeries,  the  amplitude  of  the  wave  would  be  inversely  proportional  to  the  mode 


number  cubed. 


where 


n,  m,  q nmq 


n _ m _ q 
b^a 


P a 

A » -L- 

»,m.q 


Applying  the  dynamic  fracture  criterion  to  the  amplitude  of  the  wave,  we 


obtain 


P a 

— ~ - S 

mb  * 
c 


/ P Y/3  / x 1/3 

'••M  ft)  ■ 


The  fragment  size  is  a/m  or 

c 


/2SsV/3  / 2 \1/3  (P 
ly d = hrj  (a  b)  £or  \ tt 


where  the  constant  was  chosen  to  be  continuous  with  the  one -dimensional 
results  when  i = a.  For  1 >a,  we  again  UBe  the  one  dimensional  result. 


ft) b for  ft  s 0 


Hence,  the  three-dimensional  analysis  is  appropriate  only  if  the  fragment 
size  is  less  than  a.  Otherwise,  shattering  is  one -dimensional  along  the  length 
of  the  specimen.  -23- 


HI.  ICE  CRYSTAL  BREAKUP 


St 


A)  Der.cripti.on  of  the  Model 

A ■ the  fractured  ice  crystal  passes  through  the  vehicle  shock 
layer,  it  will  break  up  and  deform  by  th«  acHor  of  the  shock  layer  gas. 

The  deformation  of  the  fragment  cloud  is  studied  in  section  IV.  However, 


on  a time  scale  small  with  respect  to  the  cloud  deformation  time,  ice 
crystal  breakup  must  be  assessed  in  order  to  determine  the  appropriate 
model  for  the  fragment  cloud.  That  is,  does  the  fragment  cloud  behave  as 
a separate  fluid,  a two  phase  flow  or  as  a fluidized  bed? 

To  begin  to  answer  this  question,  we  see  from  Section  II-B  that 
the  maximum  velocity  imparted  to  the  solid  particle  is 


a A. 
o k 


a P 

O B 


5 3 

Since  a ~ 3 km/sec,  E ~1Q  atm  and  P ~ 10  atm  (max),  we  obtain 
o s 

U ~ 30m/sec 
max 

The  expansion  Mach  number  (relative  to  the  gas)  corresponding  to  U 

max 

would  be  0.  1 upstream  of  the  vehicle  shock  and  0.  01  downstream  of  the 
shock.  Hence,  the  expansion  is  much  too  slow  to  be  thought  of  as  an 
"explosion".  Rather,  we  think  of  it  as  a "convective  separation"  where  the 
crystal  fragments  separate  with  some  characteristic  velocity  indicative  of 
the  initial  fracturing  process.  As  the  fragments  separate,  shock  layer  gas 
will  influx  through  the  cracks  as  shown  in  Figure  8.  The  time  scale  for  gas 


to  penetrate  the  fragment  cloud  is  short  compared  to  the  time  required 


mctA«ppkk(k 


I 

I 


I 


Fiq.  8 tee  Crystal  Breakup 
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for  the  fragment  to  respond  dynamically  to  the  gee.  In  fact,  the  latter 
time  scale  is  the  deformation  time  and  this  is  the  distinction  between 
crystal  breakup  and  cloud  deformation. 

As  the  gas  penetrates  through  the  cracks,  the  gas  seeks  pressure 
equilibrium  {^7P  -*0).  However,  the  pressure  difference  across  a frag- 
ment is  i^F,  which,  if  in  excess  of  will  fracture  the  fragment.  Con- 
tinuous fracturing  of  the  fragments  would  increase  the  surface /volume  of 
the  solid  and  prevent  further  influx  of  gas.  If  the  incoming  gas  achieves 
pressure  equilibrium  within  the  fragment  cloud  before  further  fracturing 
can  occur,  the  fragment  cloud  must  be  treated  as  a two  phase  flow.  How- 
ever, if  further  fracturing  prevents  the  gas  from  penetrating  the  fragment 
cloud,  the  cloud  may  be  thought  of  as  a separate  gas  which  deforms  in 
much  the  same  way  that  a liquid  drop  deforms.  In  the  following  sections, 
it  is  shown  that  the  fragment  cloud  may  be  treated  as  a separate  fluid  which 
deforms  as  a continuous  medium. 

B)  Formulation  of  the  Model 

Investigating  the  stagnation  region  oi  the  fragment  cloud,  we  assume 
that  the  velocity  of  the  fragments  is  radial  with  magnitude  increasing  linearly  in 


where  R is  the  initial  radius  of  the  fragment  cloud  or  ice  crystal.  If  the 
o 
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where  r(o)  is  the  independent  coordinate  r end  r(t>  is  the  coordinate  of 
the  fragment  that  was  initially  at  r.  The  porous  volume  within  the  frag- 
ment cloud  may  be  expressed  as 


* (inr5,,)  . ±"r3)  § 


where  Q is  the  solid  angle  aubt-^nded  by  the  stagnation  area.  For 


t <<  R /U 
o 


°*oP.*  3 

lCR-r 


dV  * A dr 
P P 


(8) 


and 

3fla  P tr 

A = — ? 

P ERo 

The  equations  of  motion  for  the  gas  are  written  using  a control  volume 
fixed  to  the  fragments.  It  is  assumed  that  the  velocity  of  the  fragments  is 
constant  in  time  and  small  with  respect  to  the  gas  velocity. 

Gas  Continuity: 


Gas  Momentum : 

6 — 6 / 2 \ bp 

(P  u dV  ) + v—  ( P u A ) dr  = -A  —■  dr  - FdA  NdA 
Ot  p 0 r V p/  p Or  s n 


where  P is  the  gas  density,  u the  gas  velocity  and  p is  the  gas  pressure. 
The  shear  force  between  the  fragments  and  the  gas  is  denoted  by  F and  the 
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normal  pressure  force  is  denoted  by  N. 


To  determine  F,  consider  the  incompressible  flow  between  parallel 
plates  separated  by  distance  2h 

n.2  2w  V2 

U = U (h  - y )/h 

h 

where  y is  measured  from  the  centerline,  equidistant  between  the  plates. 
The  average  gas  velocity  u is  related  to  u by  continuity 


u = f " ( ')  = — u 

\2h;  3 u< 


The  shear  stress  F is  determined  from 


- 6u 

F * 

ay 


,—  77 

O.  3u  H 

h ' h 


-h 


where  U is  the  gas  viscosity. 

To  determine  h(r,t),  note  that  the  p».-rous  volume  per  fragment  is 

6h/  where  -1  is  the  fragment  size.  The  number  of  fragments  in  volume 
2 2 3 

element  Or  dr  is  Or  dr/-f  . Therefore,  the  porous  volume  element  is 


dvP  = 


H)  (“H 


Comparing  to  eq  (8),  we  obtain  h. 


h = 


a Pit 

o s 

2ER 


The  surface  area  of  each  fragment  is  6<l  . Thus,  the  net  surface  area  in 


volume  element  Or  dr  is 


dA 


= 612 


(‘ r ‘) 


-29- 


FdA  = 

■ 


36  H u ER  Or  dr 
o 


i2*  p t 


The  normal  pressure  force  N is 


N = -^  P u sign  ft). 


The  normal  surface  area  associated  with  each  fragment  is  2h3,  Thus, 


4A„  = 2hi 


wo 


— — 2 _ Z 

a P t P u sign  (u  )q  r dr 

NdA  = -2— ^ 

n 2ER  3 

o 

Using  V^,  A^,  FdA^  and  NdA^  , the  continuity  and  momentum  equations 


become,  respectively. 


hh  <Pu'J)  <”>  * 0 


2 2 2 

ft-  - 12  u HE  R ~~ 

— Ou  , — — Ou  6n  o Pu  — 

P Ft  U $7  6?  <T  ~ ?,?  2 + TT  818n  * u ) = 0.  (10) 

a P * t 
o s 

The  fragment  size  1 appearing  in  eq,  (10)  is  a time  dependent  quantity  which  is 
yet  to  be  determined.  If  we  apply  a pressure  gradient  across  » crystal  frag- 
ment, that  fragment  will  fracture  when 


|*  > o 

Or  s. 


The  shattering  theory  indicates  that  the  fragment  will  break  into  smaller 

- 1/3 

fragments  of  size  (o  / l z£  ) t in  a time  scale  of  order  3/a  . 

* • Or  o 


i 


Therefore 


or 


(ID 


Equation*  (9),  (10)  and  (ll)  represent  three  equation*  in  the  four  unknowns 
P»  u,  p and  l.  In  the  next  two  sections,  they  are  solved  by  specifying  P in 
the  incompressible  and  sonic  limits. 

C)  Incompressible  Solution 

In  the  incompressible  limit,  the  continuity  equation  yields  u, 

u = - r/Jt 


and  the  momentum  equation  yields  p (*.r.t) 


sj  4U  EVr 

6r " , 2P  ^t3 

O S 


- - 2 

4 p r Pr 

9 12  5^.t2 


!12) 


Substituting  -JE  into  the  shattering  relation,  eq.  (11),  yields 
Or 


b(l/i  ) 

o __ 

6T 


l + 


W T 


AK)  . (A^f  „3,  3 (U(^) 


- K3^3 


-1/3 


(13) 
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where  1 


■ft) 


i/3 


* = t/(i  /a  ) 
o o 


Pa  R 
Rey  = ■■  _°  ° 
P 


• (*>,n(£r(*.rw 

8 


•(!)' 


/3p.2/V/3/W/3p.) 


O \ 1/9 


5/9 


K = (24) 


1/3  (*) 


For  the  ice  cyratal/ shock  layer  interaction  problem,  we  are  interested 
in  the  following  range  of  values. 

Rey  = 10*  to  10^ 

P =10  atm  to  ID3  atm 


o =10  atm 

8 • 


P = 100  P 


E = 10  atm 


Figure  9 illustrates  that  our  range  of  interest  (cross  hatched  region)  lies  in 

the  domain  where  £ >>£,,£  »K£.,  and  £ >>1.  Thus,  for  r near  R , 

o 1 o i o o 

equation  (13)  becomes 
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Rey  * p o0 Mm 


rig.  9 Ordering  Of  The  Shattering 
Equation  (incompressible  Flow) 
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6(i/i  ) 
o 

&T 


= -l  + 


ar 


€ 


o 

The  nature  of  equation  (14)  depend*  on  the  size  of  6 . For  € «l, 

o o 

chattering  will  stop  at  some  value  of  X near  f (shattering  stops  when 

ax  „ , 

— ^ * 0 or  when  * = 0).  However,  for  the  case  of  interest,  (eQ>>l) 
the  approximate  solution  to  (14)  is 


(14) 


(15) 


for  ts  i (X  = 0 for  T»  i ).  Thus,  continuous  shattering  at  r = R breaks 

o 

the  outer  fragments  into  tiny  fragments  which  will  prevent  further  penetra- 


tion by  the  shock  layer  gas. 

To  estimate  the  penetration  distance  of  the  shock  layer  gas,  we 
assume  that  the  pressure  gradient  at  r = Rq  is  continuous  over  the  penetra- 
tion distance,  6. 


6 a P f fx2)  f6  <<R  ) 0* 

s \ w r / \ max  o / 

From  (15)  and  (IE),  the  pressure  gradient  is  expressed  in  terms  of  dimen- 
sionless parameters  previously  defined. 


&P 
6 r 


3 

c c 

^ o a 

i t3(1.  T) 
o 


2 


The  penetration  distance  becomes 


(r  as  R ) 

o 


_6_ 

X 

o 


PbT3(i_T)2 

a 

C 'a 

o s 
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or 


o 30 cG 

8 O 

For  the  range  of  parameters  illustrated  in  Figure  9,  ^max  much  less 
than  ^ . Thus,  the  shock  layer  gas  does  not  penetrate  the  fragment  daid,  and 
the  fragment  cloud  may  be  treated  as  a separate  gas  or  incompressible 
fluid  without  surface  tension. 

D)  Sonic  Solution 

The  gas  flux  into  the  porous  regions  of  the  fragment  cloud  iB  more 
Hkely  to  be  sonic  than  incompressible.  Rewriting  the  shattering  equation 


for  u = u , we  obtain 


(f) 

■ = -1  + 


Mr* 


(i) 


where  we  assumed 


*_  * 

P u u = F 


s)'4 


c = 

2 


a fi  P 73 
s s 


T* 

For  u = a. ^ and  the  range  of  parameters  indicated  in  Section  C),  Figure  10 
illustrates  that  the  C_  term  in  eq.  (17)  is  dominant. 


b(1/ 1 


= - 1 + 


(r 


3t2/3 
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For  C2  >>1 


I" 

o 


=?  1 - T 


and  the  pressure  gradient  is  expressed  as 


ii  - CA 

' t t2(1_t,2 
o 


(18) 


The  penetration  depth  of  the  shock  layer  gas  is  determined  from  eqs.  (16) 
and  (18) 


b_ 

l 

o 


P8T2(i_T)2 

C 3a 

2 s 


or 


max 


16C  3o 

2 s 


For  the  range  of  parameters  indicated  in  figure  10,  6 <<t  . The  shock 

max  o 

layer  gas  does  not  penetrate  the  fragment  cloud  and  the  fragment  cloud  may 
be  treated  as  a separate  gas  or  incompressible  fluid  without  surface  tension. 


IV.  FRAGMENT  CLOUD  MOTION 


A)  Model  For  the  Cloud  Motion 

The  analysis  of  the  ice  crystal  breakup  indicates  that  we  may 
treat  the  fragment  cloud  as  a liquid  drop  without  surface  tension.  We 
treat  only  the  mean  deformation  of  the  cloud,  and  any  consideration  of  aero- 
dynamic stripping  or  catastrophic  breakup  is  omitted.  When  a liquid  drop 
is  out  of  dynamic  equilibrium  with  a gas,  the  stagnation  pressure  will  tend 
to  flatten  the  liquid  drop,  increase  its  drag  area  and  cause  it  to  accelerate 
faster  than  a rigid  particle.  This  effect  would  tend  to  slow  down  shattered 
ice  crystals  that  penetrate  the  shock  layer  of  a hypersonic  vehicle  and 
thereby  help  "shield"  the  vehicle  from  ice  erosion. 

To  model  this  effect,  we  extend  an  approach  used  by  Reinecke  and 
Waldman^  in  their  rain  drop  studies.  Writing  the  drop's  momentum 
equation  lateral  to  the  relative  gas  velocity,  Reinecke  and  Waldman^  use 
a moment  method  to  develop  an  equation  for  the  lateral  deformation.  Assum- 
ing that  the  lateral  liquid  velocity  is  linear  in  the  lateral  coordinate,  they 
obtain 


+ 2 


= 4C 


P 


* 


MWW 


where 


r 


t. 


t 


A = R/R  or  D/D 
o o 

„ vVc- 
T = — — 

o 

= initial  relative  gas  velocity 

p J 

C = gas  density/liquid  density  = /P 
C = pressure  coefficient  at  Raindrop  Stagnation  Point 

The  Reinecke  and  Waldman  approach  is  limited  to  times  for  which  the  rela- 
tive velocity  u is  approximately  equal  to  its  initial  value  u^.  Defining  the 
dimensionless  relative  velocity  v 


v 


we  extend  the  Reinecke- Waldman  equation  to  include  a time  dependent 
relative  wind 


d2A 

dT2 


and  express  the  dynamics  of  u 


- (Sr-*.. 


(19) 


/4  3\  du 

p(lnRo  ) dr  = 


1-2  2 
- 2 P U CDT7R 
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1 


or 


£ ■ - 


2 2 
v 


(20) 


A simultaneous  solution  to  eqs.  (19)  and  (20)  is  needed  to  determine  the 
slowdown  and  lateral  spreading  of  the  deformable,  incompressible  drop. 
B)  Early  and  Late  Time  Solutions 

To  obtain  a solution  eqs.  (19)  and  (20),  it  is  natural  to  seek  an 
early  time  solution  of  the  form 


VC  Vj 


V * V + 

o 


A = A A . - 

o 1 


as  £ "*  0 
as  C **0 


where 


and 


A =\2C  T + constant  (Reinecke- Waldman  Result) 
o p 


v = 1. 
o 


Obtaining  the  next  order  term  in  v. 


C T 
P 


we  note  that  the  expansion  is  singular  when  both  T and  A are  0 (£ 
Rescaling  the  problem  to  "late  time",  we  define 

£ > e 1/6t 

1 th 

r = c a 


-1/6. 
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and  the  equations  becom-s 


rdz 

d£2 

dv 

= 

The  early  time  sobucion  provides  the  initial  conditions  to  (Zl)  and  (ZZ). 

vk»  = 1 

r ( 0)  =0  as  C **0 


m.  -£T 

’ p 

Numerical  solutions  for  v^^aodr^^are  illustrated  in  Figures  11  and 
1Z  respectively.  Over  the  ;:ime  scale  a substantial  deviation  from 
v = l is  noted.  The  decrease  in  v is  enhanced  by  the  deformation  of  the 
drop.  This  deformation  causes  the  liquid  drop  to  accelerate  more  rapidly 
than  a rigid  particle. 

If  we  imagine  a ahook  tube  experiment  where  a stationary  liquid 

drop  is  suddenly  exposed  to  a relative  wind  u^,  the  drop  displacement 

(x„,J  in  the  shock  tube  becomes 
' ST 


dx 


ST 


= U - U 
7 


(23) 


In  the  "Reinecke- Waldman  Coordinate",  we  nondimensionalize  to  Dr 
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t# Wrjii'.i  fl  ,A/ «•  -<M*  >4  tlt.'fii 


Fig.  12  Late  Time  Solution  For  The 
Lateral  Deformation 


NONDIMENSIONAL  TIME  <€*  U2t/D0) 


XST  = *ST,Dc 


wfatms  in  the  present  coordinates,  X4>T  is  re- scaled  by  r 

>.2/, 

Y__  = € 3 X0_ 

ST  ST 


Eq.  (23)  becomes 


= 1 - v 


A numerical  solution  of  eq.  (24)  (coupled  to  (21)  an'4  (22)  ) is  illustrated 
in  Figure  13.  The  results  illustrate  a very  rapid  acceleration  at  small  time. 
Making  a small  VI  expansion  of  (21),  (22)  and  (24),  we  obtain 


=^jzc  H 


- IcrS^5 


VST  * I SCp£4 


Rescaling  to  the  Reinecke-  Waldman  coordinates 


l 4 

XST  "5  ^ (Liquid  Drops) 


whereas,  for  so5id  particles,  we  integrate  eq.  (20)  with  A *=•  l. 


V - 2 r T2 
ST  8 D 


(Solid  Particles) 
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SHOCK  TUBE  DISPLACEMENT  (X8T«VD0) 


Thus,  the  acceleration  of  a liquid  drop  is  significantly  different  from 

that  of  a solid  particle  and  the  difference  in  demonstrated  empirically  in 

Figure  14.  Recent  data  of  Reinecke  and  Waldman^  ^ together  with  that 
(2)  4 

of  Engel  closely  follows  the  T curve  except  for  X _ > 10, where  the 

b JL 

4 

T results  must  yield  to  the  numerical  solution.  Slightly  smaller  values 
of  Cp  or  Cp  would  yield  even  better  agreement  with  the  data. 

C)  Application  to  Stagnation  Point  Erosion 

Having  developed  a model  for  the  dynamics  of  a deformable, 
incompressible  fragment  cloud,  we  will  now  apply  the  theory  to  the  stagna- 
tion point  of  a hypersonic  vehicle  and  determine  the  conditions  tinder  which 
the  shattering  of  the  ice  crystal  into  a fragment  cloud  will  help  shield  the 
vehicle  from  ice  erosion. 

Considering  the  stagnation  point  shock  layer  to  be  a static  gas,  the 

penetration  (x  ) of  the  fragment  cloud  into  the  shock  layer  may  be  written 
sp 


dx 

= u 

dt 


dY 

— = v.  (25) 

Coupling  eqs.  (25)  to  eqs.  (21)  and  (23),  numerical  results  are  illustrated 

in  Figure  15.  Figures  11  and  15  may  be  used  to  determine  v(x#p).  At 

x = A , the  shock  stand  off  distance,  v represents  v , the  relative  velo- 
sp  s 1 

2 

city  at  impact,  and  v represents  the  relative  kinetic  energy  at  impact. 
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XST- SHOCK  TUBE  DISPLACEMENT  (XST/O0) 


rig.  14  Dynamics  Of  Liquid  Drops 


z 
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NONDIMEN  SIGNAL  TIME  ( Uzt/D0) 


To  determine  for  a solid  particle. 


use  Eq.  (20)  with  A = 1 (D  = D ). 

o 


Integration  yields 


and 


f 

t 


At  x ~ A , v = v,,  therefore 
*p  s I 

Vj  = exp  (-iC^C  Ag  M.  (Rigid  Particles) 

The  relative  impact  kinetic  energy  of  solid  particles  ie  compared 
to  that  of  a fragment  cloud  of  the  same  initial  size  in  Fig.  16.  The  results 
indicate  that  the  smallest  size  fragment  cloud  that  will  not  cause  damage 
to  a hypersonic  vehicle  may  be  an  order  of  magnitude  greater  than  the 
minimum  tolerable  rigid  particle.  A specific  case  s eated  by  Wu1^  is 
illustrated  in  Big.  17  and  compared  to  the  present  results.  Wu  com- 
pared the  results  foi  a rigid  particle  to  that  for  the  melting,  vaporizing  parti- 
cle and  showed  that  the  critical  size  for  causing  damage  war  above  10fj. 

When  ice  crystal  shattering  is  accounted  for,  the  vehicle  under  consideration 
car.  tolerate  up  to  50  p diameter  crystals.  The  -*hielding  effect  obviously 
increases  rapidly  with  increasing  nose  radius.  Consideration  of  fragment 
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cloud  melting,  vaporization,  (tripping  and  Taylor  i ratability  should 
further  increase  the  i.'inimum  tolerable  size  ice  crystal. 
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V.  CONCLUSIONS 


The  problem  of  ice  crystal  impact  with  a hypersonic  vehicle  has 
been  addressed.  The  conditions  under  which  the  ice  crystal  will  fracture 
at  the  vehicle  shock  have  been  determined.  A dynamic  fracture  criterion, 
which  is  subject  to  experimental  verification,  has  been  developed  to  deter- 
mine the  sise  and  velocity  of  the  crystal  fragments.  Using  die  initial  frag- 
ment sice  and  velocity  from  the  model  for  dynamic  fracture,  it  has  been 
shown  that  the  outer  fragments  continue  to  fracture  and  prevent  the  shock 
layer  gas  from  penetrating  the  fragment  cloud.  Thus,  the  fragment  cloud 
may  be  treated  as  a separate  incompressible  fluid  without  surface  tension. 
The  dynamics  of  the  fragment  cloud  is  then  modeled  by  extending  the 

(5) 

Reinecke- Waldman  model  for  liquid  drop  deformation.  The  lateral 
spreading  of  the  fragment  cloud  may  significantly  reduce  the  ice  crystals 
impact  velocity  on  the  hypersonic  vehicle.  It  is  shown  that  the  minimum 
tolerable  fragment  cloud  may  be  an  order  of  magnitude  greater  in  size  than 
the  minimum  tolerable  rigid  particle  impacting  the  vehicle  along  the  same 
stagnation  point  trajectory.  Thus,  the  fracturing  of  the  ice  crystals  at  the 
vehicle  shock  may  significantly  reduce  the  erosion  of  the  vehicle. 
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